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Abstract;  A  numerical  model  study  of  the  Hoosic  River  Flood-Control 
Channel  was  conducted  to  determine  the  flow  conditions  with  the  channel 
as  it  presently  exists  as  compared  with  the  as-built  conditions.  Sediment 
has  accumulated  in  certain  portions  of  the  channel  resulting  in  relatively 
large  areas  of  deposited  material  that  may  inhibit  flood-flow  conveyance. 
The  two-dimensional  depth-averaged  module  of  the  Adaptive  Hydraulics 
(ADH)  finite  element  flow  solver  was  used  to  obtain  velocity  information 
and  water-surface  elevations.  Validation  of  the  numerical  modeling  system 
was  completed  by  comparing  simulation  results  with  published  physical 
model  data.  The  model  was  then  used  to  evaluate  the  flow  conditions 
associated  with  the  as-built  channel  configuration.  The  existing  channel 
conditions  were  then  simulated  to  determine  the  water-surface  elevations 
that  are  to  be  expected  with  design  discharge  and  the  channel  in  the 
present  day  condition.  Finally,  the  flow  conditions  with  two  alternative 
channel  modifications  were  evaluated.  These  modifications  were  pro¬ 
posed  channel  restoration  configurations  designed  to  improve  habitat  and 
aesthetics. 
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graphs  in  this  report  are  taken  from  their  report. 

COL  Gary  E.  Johnston  was  Commander  and  Executive  Director  of  ERDC. 
Dr.  James  R.  Houston  was  Director. 


ERDC/CHLTR-10-1 


ix 


Unit  Conversion  Factors 


Multiply 

By 

To  Obtain 

cubic  feet 

0.02831685 

cubic  meters 

feet 

0.3048 

meters 

miles  (U.S.  statute) 

1,609.347 

meters 

ERDC/CHLTR-10-1 


1 


1  Introduction 

Background 

The  Hoosic  River  begins  in  western  Massachusetts  and  flows  generally 
north  and  northwest,  crossing  the  southwest  corner  of  Vermont  to  the 
vicinity  of  Eagle  Bridge,  New  York,  where  it  turns  west  and  joins  the 
Hudson  River  about  i6  miles  above  Troy,  New  York  (Figure  i).  It  is 
66  miles  long  and  has  a  drainage  area  of  713  square  miles.  Adams,  Massa¬ 
chusetts  is  situated  on  the  south  branch  of  the  Hoosic  River,  about  8  miles 
below  the  source  of  the  river.  The  flood-control  project  for  Adams  provides 
a  man-made  channel  through  the  town.  The  upstream  reach  of  the  project 
consists  of  a  rectangular  concrete  channel  varying  from  40  to  75  ft  in  width 
(Figure  2).  The  downstream  reach  is  a  trapezoidal  earth  channel  having  a 
base  width  varying  from  65  to  90  ft. 

Purpose 

The  purpose  of  this  study  was  to  determine  the  flow  conditions  in  the  tran¬ 
sitions  and  curves  in  the  supercritical  channel,  at  the  confluences  with  the 
Fiske  Brook  and  Tophet  Brook  Channels,  at  the  drop  structures,  and  in  the 
reach  where  the  flow  transitions  from  supercritical  to  subcritical  flow.  Par¬ 
ticular  attention  was  directed  to  evaluating  the  consequences  of  the  chan¬ 
nel  changes  that  have  occurred  since  initial  construction.  These  changes 
are  predominately  the  deposition  of  sediment  in  the  lower  portion  of  the 
modeled  reach.  The  study  was  primarily  for  the  purpose  of  evaluating  vari¬ 
ous  modifications  to  the  lower  end  as  part  of  a  channel  restoration  project. 
The  model  was  to  determine  if  restoration  alterations  would  change  water- 
surface  elevations  associated  with  the  design  discharge. 

Approach 

The  first  phase  of  the  study  tested  the  numerical  model’s  ability  to  repro¬ 
duce  the  hydraulic  conditions  observed  in  a  previously  published  physical 
model  study.  After  the  numerical  model  was  validated,  the  as-built  condi¬ 
tions  were  simulated  to  establish  water-surface  elevations  associated  with 
the  initial  design  and  construction.  After  as-built  flow  conditions  were 
established,  flow  conditions  produced  with  the  existing  conditions  of  the 
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Figure  1.  Location  of  model  studies. 
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Figure  2.  Hoosic  River  through  tow/n  of  Adams,  MA. 


channel  were  computed.  The  existing  conditions  include  large  deposits  of 
sediments  at  various  locations  within  the  channel  (see  Figures  3  and  4),  so 
there  was  a  question  as  to  how  this  may  change  water-surface  elevations 
during  flood  events.  The  final  phase  of  the  study  was  the  hydraulic  evalu¬ 
ation  of  restoration  alternatives  to  determine  their  impacts  on  water- 
surface  elevations. 
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Figure  3.  Sediment  deposition  and  vegetation  at  confluence  of  Hoosic  River 

and  Tophet  Brook. 
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2  Validation  Modei 

Background 

The  two-dimensional  (2-D)  model  was  tested  to  see  how  well  it  could 
reproduce  the  water-surface  elevations  within  a  concrete  high-velocity 
channel  having  complex  geometric  features  including  spiral  curves,  invert 
banking,  and  very  few  tangents. 

Validation  required  a  data  set  of  observed  flow  conditions.  A  portion  of  the 
Hoosic  River  flood-control  project  through  North  Adams,  MA  (Figure  i) 
was  previously  tested  in  a  physical  model  (U.S.  Army  Engineer  Waterways 
Experiment  Station  1952).  The  2-D  numerical  model  was  validated  by 
comparing  computed  water-surface  elevations  with  those  measured  in  the 
physical  model.  The  actual  portion  of  the  Hoosic  River  that  was  modeled 
for  validation  purposes  is  just  as  complicated;  perhaps  more  so,  than  the 
study  reach  through  Adams.  The  project  was  a  good  case  to  test  the 
numerical  model’s  ability  to  simulate  the  complicated  flow  throughout  the 
study  reach.  The  validation  (North  Adams)  reach  contained  23  curves  with 
compound  spirals  and  invert  banking,  a  confluence  of  two  supercritical 
channels,  and  a  twin  stilling  basin  used  to  contain  hydraulic  jumps  prior  to 
introducing  the  flow  to  a  soft-bottom  channel.  The  longest  tangent 
between  compound  curves  was  less  than  5  times  the  channel  width.  Most 
curves  are  adjoined  in  that  the  spiral-to-tangent  (ST)  ending  one  curve  and 
the  tangent-to-spiral  (TS)  beginning  the  next  curve  share  the  same  station. 

The  validation  project  is  the  North  Branch  and  South  Branch  within  North 
Adams  (see  Eigure  5).  Both  channels  have  a  channel  width  of  45  ft  along 
most  of  their  length.  The  South  Branch  model  was  more  that  4,000  ft  long 
and  fell  more  than  20  ft  in  elevation.  The  slope  varied,  but  the  overall 
slope  of  the  South  Branch  was  0.005.  The  average  Eroude  number  for  this 
reach  was  1.7  and  at  the  discharge  of  10,300  cfs,  the  normal  depth  was 
8.4  ft  and  the  cross-sectional  averaged  velocity  was  27  fps.  Normal  depth 
rarely  occurs  along  the  reach  of  the  channel  due  to  numerous  curves,  but 
knowing  that  the  velocity  was  expected  to  be  about  27  fps  at  a  depth  of 
8.4  ft  provided  a  general  description  of  the  flow  conditions. 
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Figure  5.  Model  validation,  physical  and  numerical  model  limits. 
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The  North  Branch  dropped  more  than  90  ft  in  elevation  over  the  4,600  ft 
modeled,  resulting  in  an  overall  bed  slope  of  0.014.  This  very  steep  slope 
on  the  North  Branch  resulted  in  an  average  Fronde  number  of  2.8.  Con¬ 
veying  the  design  discharge  (at  the  time  of  the  physical  model  study)  of 
14,600  cfs  resulted  in  a  normal  depth  of  7.5  ft  and  a  cross-sectional 
averaged  velocity  of  43  fps  at  normal  depth.  The  actual  flow  was  always 
varying  due  to  the  changes  in  alignments,  but  knowing  that  the  velocity 
should  be  greater  than  40  fps  at  a  depth  of  7.5  ft  gave  insight  into  the  flow 
characteristics  to  be  expected. 

Physical  model 

The  physical  model  was  a  i:30-scale  flume  of  the  North  and  South 
Branches  of  the  Hoosic  River  running  through  North  Adams,  Massachu¬ 
setts.  The  model  study  was  conducted  in  increments  between  1950  and 
1957  at  the  U.S.  Army  Engineer  Waterways  Experiment  Station.  Channel 
configurations  (plans  and  profiles)  were  taken  from  the  physical  model 
report  (U.S.  Army  Engineer  Waterways  Experiment  Station  1952).  Elow 
depths  were  scaled  from  the  values  plotted  in  the  report.  The  recom¬ 
mended  design  was  reproduced  in  the  numerical  model. 

The  i:30-scale  hydraulic  model,  which  was  constructed  of  plywood,  was 
hydraulically  too  rough.  The  friction  differences  are  attributed  to  the 
differences  in  the  model  and  prototype  Reynolds  numbers.  In  an  attempt 
to  better  simulate  the  prototype  roughness,  the  slope  of  the  physical  model 
was  increased.  The  physical -model  report  mentioned  that  the  model’s  bed 
slope  was  increased,  but  it  does  not  provide  the  actual  resulting  slopes. 

Data  consisted  of  velocity  and  depth  measurements  at  various  locations 
along  the  North  and  South  Branches.  These  depth  measurements  were 
used  to  develop  water-surface  profiles  along  the  walls.  The  turbulence  in 
the  high-velocity  flow  produces  only  a  quasi-steady  flow  condition.  The 
depth  at  a  particular  location  varies  in  time.  Generally,  the  physical  model 
study  reported  maximum  values  of  depth,  although  particular  tests 
recorded  temporally  averaged  values.  Because  the  physical  model  test 
reported  maximum  depths,  it  was  expected  that  the  computed  depths, 
which  are  time  averaged  over  turbulence  scales,  would  be  somewhat  less 
than  those  recorded  in  the  physical  model  study. 
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Numerical  model 

The  channel  configuration  is  shown  in  Figure  5.  More  than  6,000  ft  of  the 
North  Branch  was  reproduced  from  sta  205+96.03  to  the  confluence  with 
the  South  Branch  at  sta  141+75.  Forty  two  hundred  feet  of  the  South 
Branch  was  reproduced  from  the  upstream  limit  at  sta  183+75  to  the 
confluence  at  sta  141+75.  The  main  channel  reach  extended  from  down¬ 
stream  of  the  confluence  to  the  stilling  basin,  a  total  of  1,210  ft. 


The  shallow-water  (or  long-wave)  equations  are  a  result  of  the  vertical 
integration  of  the  equations  of  mass  and  momentum  conservation  for 
incompressible  flow  under  the  hydrostatic  pressure  assumption.  The  flow 
depth  (/i),  the  x-component  of  velocity  (u),  and  the  y-component  of 
velocity  (u)  define  the  dependent  variables  of  the  fluid  motion.  If  the  fluid 
pressure  at  the  surface  is  taken  as  zero  and  the  free-surface  stresses  are 
neglected,  the  shallow- water  equations  are  given  as: 


dU  dF  dG 

- 1 - 1 - h 

df  dx  dy 


H  =  o 


where: 


h 


U 


uh 


vh 


(1) 


(2) 


F 


hu 

hu^  +^gh^ 

P 

huv  —  h—^ 


(3) 


P 
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hv 


G 


huv  —  h^^ 

P 

hv^+^gh^-h^ 

P 


(4) 


and: 


H 


0 


dz^ 


n^g 


u^u^+v^ 


dy 


2  vgu 
^jfg 


(5) 


here: 


p  =  fluid  density 
g  =  gravitational  acceleration 
Zb  =  channel  bed  elevation 
n  =  Manning's  roughness  coefficient 

Co  =  a  dimensional  constant  (Co  =  i  for  SI  units  and  2.208  for 
U.S.  Customary  units) 

and: 

o 's  =  Reynolds  stresses  due  to  turbulence,  where  the  first  subscript 
indicates  the  direction,  and  the  second  indicates  the  face  on 
which  the  stress  acts. 

The  equations  are  discretized  using  the  finite  element  method  in  which  u, 
V,  and  h  are  represented  as  linear  polynomials  on  each  triangular  element. 
The  numerical  model  is  the  Adaptive  Hydraulics  (ADH)  code 
(http://adh.usace.armv.mil).  This  code  features  mcsh  adaption  whereby  the  com¬ 
putational  mesh  is  automatically  refined  in  areas  where  needed  to  provide 
an  accurate  solution. 
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The  computational  mesh  developed  to  model  the  validation  flow  condi¬ 
tions,  consisted  of  4548  nodes  and  7963  triangular  elements  (Figure  6). 


The  majority  of  the  mesh  was  built  with  elements  that  were  approximately 
33  ft  by  11  ft  in  the  longitudinal  and  lateral  direction,  respectively.  Greater 
resolution  (smaller  elements)  was  used  in  the  vicinity  of  the  submerged 
weir  in  the  northern  branch,  the  confluence  of  the  northern  and  southern 
branches,  and  the  spillways  in  the  main  channel.  ADH  was  set  to  refine  the 
mesh  one  level.  Reflnement  automatically  splits  the  elements  where  the 
computed  residual  exceeds  a  specifled  threshold. 

The  tailwater  was  set  at  el  677. at  the  outflow  boundary  which  was  at  the 
western  end  of  the  main  channel.  The  flow  at  the  upstream  end  of  each 
branch  was  subcritical  so  a  subcritical  inflow  boundary  condition  was 
specified  for  each  branch.  Inflow  discharges  of  14,600  cfs  and  10,300  cfs 
were  applied  to  the  boundaries  of  the  northern  and  southern  branches, 
respectively.  The  flow  accelerated  to  supercritical  due  to  the  steepness  of 
the  channel  slope.  The  transition  from  subcritical  to  supercritical  flow 
occurred  in  the  channel  contraction  on  the  South  Branch  and  at  the  sub¬ 
merged  weir  on  the  North  Branch. 


1  All  elevations  (el)  are  In  feet  referenced  to  the  National  Geodetic  Vertical  Datum  of  1929. 
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The  model  was  run  from  an  initial  condition  assuming  zero  velocity  and 
uniform  depth  to  steady  state  by  advancing  in  time.  A  Manning’s  coeffi¬ 
cient,  n,  was  set  in  the  model  based  on  channel  conditions.  Each  channel 
was  given  an  n  value  of  0.016. 

The  water-surface  profile  was  obtained  from  the  numerical  model  results 
by  extracting  depth  values  from  the  steady-state  solution  at  the  boundary 
nodes  along  the  model  walls.  These  values  were  then  plotted  and  com¬ 
pared  to  the  profiles  published  in  the  physical  model  report  (U.S.  Army 
Engineer  Waterways  Experiment  Station  1952). 

Validation  results 

The  primary  objective  of  a  flood-control  channel  is  containing  flood  flows. 
So,  the  primary  objective  in  modeling  these  channels  was  accurately  com¬ 
puting  the  depth  of  flow  throughout  the  flow  domain.  Comparisons  of  the 
physical  model  data  with  numerical  model  results  provided  an  evaluation 
of  the  numerical  model’s  ability  to  reproduce  the  observed  flow  depths. 

Physical  model  invert  elevations  along  the  right  and  left  bank  were  taken 
from  the  physical  model  report.  Depths  were  scaled  from  the  report  and 
transferred  to  water-surface  elevations  using  the  invert  elevations.  Both 
water-surface  and  bed  elevation  profiles  were  then  plotted  at  an 
exaggerated  (10:1)  vertical  scale. 

The  validation  run  was  performed  to  compare  the  numerical  computations 
to  the  flow  depths  measured  in  the  physical  model.  Examples  of  results 
from  the  numerical  model  are  provided  in  Eigures  7-9.  Comparison  of  the 
numerical  model  results  and  the  physical  model  data  are  provided  in 
Eigures  10-23.  These  figures  show  water-surface  profiles  along  the  phys¬ 
ically  modeled  reaches  of  the  North  and  South  Branches  of  the  Hoosic 
Rivers  in  North  Adams.  Each  figure  provides  both  observed  and  computed 
profiles  along  the  right  and  left  walls  (referenced  to  looking  downstream). 

The  figures  illustrate  that  the  model  accurately  reproduced  the  water- 
surface  elevations  along  the  walls  of  the  North  Branch.  The  computed 
depths  are  slightly  lower  than  those  recorded  in  the  physical  model  along 
the  lower  end  of  the  South  Branch  (sta  169+00  to  148+00).  However,  the 
differences  are  small  and,  as  previously  noted,  the  physical  model  tests 
recorded  maximum  depths  in  the  turbulent  supercritical  flow.  The  model 
does  simulate  the  super-elevation  of  the  water  surface  in  the  banked  bends. 
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These  simulation  results  show  that  the  2-D  ADH  model  is  an  appropriate 
tool  to  use  in  evaluating  flow  conditions  in  the  Hoosic  River  Flood-Control 
Channel. 


Figure  7.  Model  validation,  numerical  model  results,  water-surface  contours  downstream 
of  the  North  and  South  Branch  confluence. 


Figure  8.  Model  validation,  numerical  model  results  at  the  confluence  of  the  North 
and  South  Branches,  looking  downstream. 
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Figure  9.  Model  validation,  numerical  model  results  at  the  confluence  of  the  North 
and  South  Branches,  looking  upstream. 
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-  WSE,  LeftV\^ll,  2D  Computational  Model  North  Branch 

NOTE:  SIDESOFCHANNELAREREFERENCEDTOLOOKING  DOWNSTREAM  STA  205+96.03  TO  193+41.59 


Figure  10.  Model  validation,  North  Branch,  plan  and  profiles,  sta  205+96.03  to  193+41.59,  left  wall. 
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Figure  11.  Model  validation,  North  Branch,  plan  and  profiles,  sta  205+96.03  to  193+41.59,  right  wall. 
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Figure  12.  Model  validation,  North  Branch,  plan  and  profiles,  sta  193+41.59  to  178+67.97,  left  wall. 
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Figure  13.  Model  validation,  North  Branch,  plan  and  profiles,  sta  193+41.59  to  178+67.97,  right  wall. 
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Figure  14.  Model  validation,  North  Branch,  plan  and  profiles,  sta  178+67.97  to  165+49.51,  left  wall. 
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Figure  15.  Model  validation,  North  Branch,  plan  and  profiles,  sta  178+67.97  to  165+49.51,  right  wall. 
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Figure  16.  Model  validation,  North  Branch,  plan  and  profiles,  sta  165+49.51  to  155+98.72,  left  wall. 
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Figure  17.  Model  validation,  North  Branch,  plan  and  profiles,  sta  165+49.51  to  155+98.72,  right  wall. 
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Figure  18.  Model  validation,  South  Branch,  plan  and  profiles,  sta  183+75.00  to  169+61.24,  left  wall. 
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Figure  19.  Model  validation,  South  Branch,  plan  and  profiles,  sta  183+75.00  to  169+61.24,  right  wall. 
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Figure  20.  Model  validation,  South  Branch,  plan  and  profiles,  sta  169+61.24  to  156+75.91,  left  wall. 
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Figure  21.  Model  validation,  South  Branch,  plan  and  profiles,  sta  169+61.24  to  156+75.91,  right  wall. 
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Figure  22.  Model  validation,  South  Branch,  plan  and  profiles,  sta  156+75.91  to  141+75.00,  left  wall. 
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Figure  23.  Model  validation,  South  Branch,  plan  and  profiles,  156+75.91  to  141+75.00,  rightwaii. 
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3  As-Built  Conditions  IVIodel 

Channel  features  and  geometry 

The  model  reproduced  the  main  channel  of  the  Hoosic  River  from  below 
Plunkett  Dam  (sta  121+00)  to  100  ft  downstream  of  Lime  Street 
(sta  0+00),  a  distance  of  12,100  ft.  The  lower  125  ft  of  Fiske  Brook  was 
reproduced.  The  Tophet  Brook  modeled  reach  extended  from  the  conflu¬ 
ence  up  to  sta  3+89.  The  detailed  junctions  where  these  channels  merge 
with  the  main  channel  were  also  included  in  the  model. 

Channel  features  include  the  weir  at  sta  113+34  where  the  flow  accelerates 
from  subcritical  to  supercritical,  three  drop  structures,  each  incorporating 
a  stilling  basin,  nine  channel  bends  that  have  spiral  curves  that  transition 
the  flow  into  and  out  of  circular  curves,  with  invert  banking  in  which  the 
invert  is  rotated  about  the  center  line.  The  spiral  curve  layout  is  in  accor¬ 
dance  with  the  Hydraulic  Design  Criteria  (Hydraulic  Design  Chart  660- 
2/2  and  sheets  660-2  to  660-2/4)  from  research  conducted  by  the  Los 
Angeles  District  (U.S.  Army  Engineer  District,  Los  Angeles,  1948).  Invert 
banking  design  was  developed  in  the  original  Hoosic  River  physical  model 
study  (U.S.  Army  Engineer  Waterways  Experiment  Station  1952).  The 
invert  banking  is  given  by: 


Az 


gr 


(6) 


here: 


Az  =  change  in  invert  elevation  at  the  wall  from  center-line  invert 
elevation 
C  =  a  coefficient 

V  =  mean  channel  velocity  of  flow  approaching  the  curve 
W  =  channel  width 
g  =  gravitational  acceleration 
r  =  radius  of  channel  center-line  curvature. 

The  coefficient,  C,  used  in  the  design  of  the  Hoosic  River  channel  in  the 
reach  through  Adams  was  1.2.  This  value  was  developed  during  the  North 
Adams  physical  model  study  (U.S.  Army  Engineer  Waterways  Experiment 
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Station  1952).  These  criteria  are  now  the  accepted  hydraulic  design  criteria 
found  in  EM  1110-2-1601. 

The  lower  reach  of  the  modeled  Hoosic  River  channel,  downstream  of  sta 
35+47.83,  is  composed  of  various  cross-sectional  shapes:  trapezoidal, 
rectangular,  and  compound  with  a  vertical  wall  on  one  side  and  a  sloping 
side  on  the  other.  This  reach  has  a  base  width  of  75  ft  and  the  side  slopes 
are  laid  iV  on  3H.  The  flow  in  this  lower  reach  is  subcritical.  The  bed  is 
composed  of  gravel,  riprap,  and  derrick  stone.  The  bends  are  simple  circu¬ 
lar  curves.  The  mild  longitudinal  slope  needed  to  sustain  subcritical  flow  is 
established  and  maintained  with  the  drop  structures  at  sta  51+50, 
sta  36+54,  and  sta  20+67. 

Computational  mesh 

A  mesh  was  developed  using  as-built  plans  (line  drawings)  provided  by  the 
U.S.  Army  Engineer  District,  New  York  (NAN).  The  line  drawings  were 
used  to  create  a  CAD  figure  of  the  channel  alignment.  The  mesh  was 
generated  from  this  CAD  drawing.  The  finite  element  mesh  had 
12,868  nodes  and  23,272  triangular  elements. 

The  computational  mesh  of  the  as-built  conditions  model  is  shown  in 
Eigure  24.  The  flow  direction  is  from  lower  left  to  upper  right  of  the  mesh. 
Details  of  the  mesh  at  each  drop  structure  are  shown  on  the  figure.  A  mini¬ 
mum  of  5  elements  were  used  to  describe  the  high-velocity  channel.  Typi¬ 
cal  element  sizes  were  approximately  15  ft  by  25  ft,  with  much  smaller 
elements  (greater  resolution)  at  junctions,  piers,  drop  structures,  and 
stilling  basins. 

Boundary  conditions  and  model  parameters 

Discharges  associated  with  the  100-year  frequency  flow  event  were  used  as 
the  design  flows.  Discharge  values  were  furnished  by  the  NAN.  Peak  flows 
for  each  channel  reach  are  shown  in  Table  1.  It  was  assumed  that  the  peak 
discharges  at  the  upstream  ends  of  the  Hoosic  River  channel,  Eiske  Brook, 
and  Tophet  Brook  occur  simultaneously,  although  the  hydrographs  show 
slight  differences  at  the  times  of  peak  at  these  inflows. 

The  inflow  boundary  condition  for  the  main  Hoosic  River  channel  was 
specified  as  subcritical.  A  discharge  of  4,271  cfs  was  applied  at  the 
upstream  end  of  the  model.  Supercritical  inflow  of  958  cfs  was  specified 
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Details  at  drop  staicture  3  and 
stilling  basin 


Details  at  drop  structure  2  and 
stilling  basin 


Figure  24.  As-built  conditions,  connputational  mesh  of  Hoosic  River  channel  and  tributaries. 
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Table  1.  Design  discharges  (100-year  frequency). 


Inflow  Boundary 

Inflow  Discharge 
cfs 

Combined  Discharge 
cfs 

Hoosic  River 

4271 

4271 

Fiske  Brook 

958 

5229 

Tophet  Brook 

1544 

6773 

at  the  upper  end  of  Fiske  Brook.  The  Tophet  Brook  inflow  was  super¬ 
critical  and  the  discharge  was  1,544  cfs.  A  subcritical  outflow  boundary 
condition  was  set  at  the  downstream  end  of  the  channel.  The  outflow 
discharge  of  6,773  cfs  was  the  total  of  the  three  inflows.  Normal  depth, 
which  is  8.2  ft  at  6,773  cfs,  was  assumed  to  occur  at  the  outflow  section 
and  since  the  centerline  bed  elevation  was  734.8  a  tailwater  of  el  743.0  was 
specified  as  the  outflow  boundary. 

The  discharges  used  to  evaluate  the  Hoosic  River  channel  were  4,271  cfs  at 
the  upstream  end  of  the  Hoosic  River  at  sta  121-1-00  down  to  the  conflu¬ 
ence  of  Fiske  Brook  near  sta  103-150.  Downstream  of  Fiske  Brook  the  flow 
rate  was  5,229  cfs  until  the  addition  of  Tophet  Brook  flow.  The  Tophet 
Brook  confluence  was  near  sta  46-150  where  the  flow  discharge  increased 
to  the  total  of  6,773  cfs. 

The  upper  reach  of  the  channel  was  concrete  and  the  lower  end  varied 
from  small  cobbles  and  derrick  stone  to  riprap.  The  bed  roughness  was 
uniform  on  the  invert  of  the  high-velocity  channel,  but  varied  along  the 
model’s  lower  reach.  A  Manning’s  roughness  coefficient,  n,  of  0.014  was 
used  for  concrete,  0.025  for  soil  and  derrick  stone,  and  0.030  for  riprap. 

As-built  conditions  results 

Flow  features  that  would  produce  adverse  flow  conditions  were  modified 
during  the  physical  model  study  reported  in  U.S.  Army  Engineer  Water¬ 
ways  Experiment  Station  (1957).  The  physical  model  study’s  recommenda¬ 
tions  were  primarily  incorporated  into  the  constructed  (as-built)  channel. 
The  numerical  model  simulation  results  confirmed  the  acceptable  perfor¬ 
mance  of  the  channel  as  originally  designed  and  constructed.  The  velocity 
in  the  high-velocity  channel  reach  varied  from  about  23  fps  to  about 
29  fps.  The  velocity  of  the  flow  over  the  drop  structure  near  sta  51+50 
was  greater  than  30  fps.  Velocities  in  the  area  considered  for  channel 


ERDC/CHLTR-10-1 


32 


modifications  are  shown  in  Figures  25-29.  Note  that  the  contour  interval 
is  different  on  each  figure.  These  contours  of  velocity  magnitude  show  the 
transition  from  supercritical  flow  upstream  of  the  first  drop  structure  at 
sta  51+50,  to  subcritical  flow  along  the  lower  reach.  The  drop  structures 
downstream  of  Tophet  Brook  maintained  the  mild  slopes  needed  to  ensure 
subcritical  flow. 

The  as-built  conditions  water-surface  profiles  are  plotted  with  the  results 
of  the  existing  conditions  which  will  be  discussed  in  the  next  section.  The 
profiles  will  be  used  to  examine  the  effects  of  the  material  that  is  deposited 
in  the  existing  channel. 
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Figure  25.  As-built  conditions,  veiocity  magnitude  contours  at  Tophet  Brook  confluence, 

sta  53+50  to  41+94.07. 
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Figure  26.  As-built  conditions,  velocity  magnitude  contours,  sta  41+17  to  29+11. 
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Figure  27.  As-built  conditions,  velocity  magnitude  contours,  sta  29+11  to  17+50. 
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Figure  28.  As-built  conditions,  velocity  magnitude  contours,  sta  20+67  to  8+38.2. 
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Figure  29.  As-built  conditions,  veiocity  magnitude  contours,  sta  9+00  to  -1+00. 
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4  Existing  Conditions  Modei 

Geometry  and  computational  mesh 

The  channel  has  changed  over  the  years  since  construction  was  completed 
in  the  1950’s.  Gravel  and  rocks  have  washed  into  the  stream  and  accumu¬ 
lated  on  the  bed.  Also,  vegetation  has  grown  up  in  the  area  of  interest.  The 
bed  elevations  in  the  as-built  conditions  model  were  changed  to  simulate 
present  day  conditions. 

The  NAN  conducted  an  engineering  survey  on  18-22  June  2007.  The  sur¬ 
vey  data  described  the  existing  bathymetry  in  the  Hoosic  River  from 
sta  1+12  to  sta  50+15.93  and  that  of  Tophet  Brook  from  sta  0+77.26  to 
sta  2+56.68.  The  data  included  bed  elevations,  photographs,  and  descrip¬ 
tions  of  the  surveyed  reaches.  The  survey  data  were  used  to  update  the  as- 
built  conditions  computational  mesh  to  reflect  the  changes  reported  in  the 
survey.  The  photographs  were  used  to  visually  check  the  mesh  for  areas 
that  had  not  changed  from  the  as-built  conditions  (such  as  concrete  chan¬ 
nel  with  no  sediment  buildup  and  drop  structure  crests).  The  pier  at  the 
Lime  Street  Bridge  has  changed  since  the  Hoosic  Channel  was  first  con¬ 
structed.  Initially,  the  bridge  had  a  single  pier  within  the  channel.  The 
bridge  was  modified  in  1996  to  free  spanning.  Therefore,  the  as-built  con¬ 
ditions  model  had  a  bridge  pier  at  Lime  Street  and  the  existing  conditions 
model  did  not. 

The  existing  conditions  mesh  had  12,255  nodes  and  22,105  triangular  ele¬ 
ments.  Element  sizes  and  mesh  resolution  were  similar  to  those  used  to 
model  the  as-built  conditions.  Details  of  the  existing  conditions  mesh  are 
shown  in  Figure  30. 

Boundary  conditions  and  model  parameters 

The  boundary  conditions  (inflow  discharges  and  tailwater  elevation)  were 
identical  to  those  used  for  the  as-built  conditions  model.  The  roughness 
coefficients  were  also  similar  to  the  as-built  conditions  model.  The  primary 
difference  in  the  setup  of  the  as-built  conditions  and  existing  conditions 
models  was  inclusion  of  material  that  had  been  deposited  within  the 
channel. 
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Details  at  drop  structure  2  and 
stilling  basin 


Figure  30.  Existing  conditions,  computationai  mesh  of  Hoosic  River  channei  and  tributaries. 
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Existing  conditions  resuits 

The  numerical  model  results  are  illustrated  on  the  three-dimensional  per¬ 
spectives  provided  in  Figures  31  and  32.  Figure  31  shows  the  super¬ 
elevation  of  the  water  surface  as  flow  passed  around  a  bend.  The  hydraulic 
jumps  downstream  of  the  drop  structure  in  the  main  channel  (near 
sta  51+50)  and  Tophet  Brook  are  shown  in  Figure  32.  The  jumps  were 
contained  within  the  stilling  basins  as  they  were  designed  to  perform. 

Velocities  in  the  lower  end  of  the  modeled  reach  where  restoration  efforts 
will  be  directed  are  shown  in  Figures  33-37.  Note  that  the  contour  inter¬ 
vals  are  different  on  each  figure.  The  range  and  intervals  of  contours  are 
different  to  accentuate  the  differences  across  the  reach  depicted  on  each 
figure. 

Detailed  water-surface  profiles  resulting  from  the  design  discharge  are 
provided  in  Figures  38-55.  The  plan  and  profile  drawings  show  the  chan¬ 
nel  alignment,  invert  elevation,  and  water-surface  profiles  for  the  as-built 
and  existing  conditions.  The  profiles  are  provided  along  both  the  right  and 
left  walls  (looking  downstream).  The  water-surface  profiles  down  to  the 
drop  structure  near  51+50  were  similar  for  the  as-built  and  existing  condi¬ 
tions.  This  is  because  the  models  were  identical  throughout  this  reach. 

This  recalculation  of  as-built  conditions  was  made  simply  for  convenience 
of  model  setup. 

The  material  deposited  on  the  bed  is  shown  in  the  bed  profiles  down¬ 
stream  of  sta  51+00.  The  existing  conditions  increase  the  water-surface 
elevation  more  than  2  ft  from  near  sta  51+50  (Figures  50  and  51)  down  to 
sta  26+00  (Figures  52  and  53).  The  water-surface  elevation  downstream  of 
sta  26+00  with  the  existing  conditions  is  more  than  1  ft  higher  than  the 
as-built  conditions  as  shown  in  Figures  52-55. 
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Figure  31.  Existing  conditions,  super-eievation  in  computed  water  surface  through  curve. 


Figure  32.  Existing  conditions,  water  surface  at  the  confluence  of  Tophet  Brook 

and  Floosie  River. 
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Figure  33.  Existing  conditions,  veiocity  magnitude  contours  at  Tophet  Brook  confluence, 

sta  53+50  to  41+94.07. 
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Figure  34.  Existing  conditions,  veiocity  magnitude  contours,  sta  41+17  to  29+11. 
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Figure  35.  Existing  conditions,  veiocity  magnitude  contours,  sta  29+11  to  17+50. 
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Figure  36.  Existing  conditions,  veiocity  magnitude  contours,  sta  20+67  to  8+38.2. 
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Figure  37.  Existing  conditions,  velocity  magnitude  contours,  sta  9+00  to  -1+00. 
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Figure  38.  Existing  and  as-built  conditions,  plan  and  profiles,  sta  121+00  to  106+73.04,  left  wall. 
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Figure  39.  Existing  and  as-buiit  conditions,  pian  and  profiies,  sta  121+00  to  106+73.04,  right  waii. 
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Figure  40.  Existing  and  as-buiit  conditions,  pian  and  profiies,  sta  106+73.04  to  93+37.43,  ieft  waii. 
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Figure  41.  Existing  and  as-buiit  conditions,  pian  and  profiies,  sta  106+73.04  to  93+37.43,  right  waii. 


ui 

o 


ERDC/CHLTR-10-1 


LEGEND 

AS- BUILT  BED 

- EXISTING  BED 

- AS-BUILT  WATER  SURFACE 

- EXISTING  WTER  SURFACE 


DISCHARGE  STATIONS 
4271  CFS  121 -too  TO  103-160 

5229  CFS  103-150  TO  46-t50 

6773  CFS  46-150  TO -1-H30 


NOTE:  SIDES  OF  CHANNEL  ARE  REFERENCED  TO  LOOKING  DOWNSTREAM 


Plan  and  Left  Side  Profiles 

STA  93-f37  43  TO  STA  79-t-40.90 


Ul 

!-» 


Figure  42.  Existing  and  as-built  conditions,  plan  and  profiles,  sta  93+37.43  to  79+40.90,  left  wall. 
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Figure  43.  Existing  and  as-buiit  conditions,  pian  and  profiies,  sta  93+37.43  to  79+40.90,  right  wali. 
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Figure  44.  Existing  and  as-built  conditions,  plan  and  profiles,  sta  79+40.90  68+13.34,  left  wall. 
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Figure  45.  Existing  and  as-built  conditions,  plan  and  profiles,  sta  79+40.90  68+13.34,  left  wall. 
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Figure  46.  Existing  and  as-buiit  conditions,  pian  and  profiies,  sta  68+13.34  to  55+25.00,  ieft  waii. 
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Figure  47.  Existing  and  as-buiit  conditions,  pian  and  profiies,  sta  68+13.34  to  55+25.00,  right  waii. 
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Figure  48.  Existing  and  as-buiit  conditions,  pian  and  profiies,  sta  55+25.00  to  41+94.07,  ieft  waii. 
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Figure  49.  Existing  and  as-buiit  conditions,  pian  and  profiies,  sta  55+25.00  to  41+94.07,  right  waii. 
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Figure  50.  Existing  and  as-buiit  conditions,  pian  and  profiies,  sta  41+94.07  to  27+19.78,  ieft  waii. 
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Figure  51.  Existing  and  as-buiit  conditions,  pian  and  profiies,  sta  41+94.07  to  27+19.78,  right  waii. 
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Figure  52.  Existing  and  as-buiit  conditions,  pian  and  profiies,  sta  27+19.78  to  12+14.29,  ieft  waii. 
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Figure  53.  Existing  and  as-buiit  conditions,  pian  and  profiies,  sta  27+19.78  to  12+14.29,  right  waii. 
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Figure  54.  Existing  and  as-buiit  conditions,  pian  and  profiies,  sta  12+14.29  to  0+00,  ieft  waii. 
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Figure  55.  Existing  and  as-buiit  conditions,  pian  and  profiies,  sta  12+14.29  to  0+00,  right  waii. 


0) 


ERDC/CHLTR-10-1 


ERDC/CHLTR-10-1 


65 


5  Restoration  Models 

Chapter  i  explained  that  the  upstream  reach  of  the  project  consists  of  a 
rectangular  concrete  channel,  whereas  the  downstream  reach  is  a  trape¬ 
zoidal  channel  with  riprap  protection.  The  trapezoidal  reach  is  a  man¬ 
made  earth  channel  that  has  a  base  width  that  varies  from  65  to  90  ft.  The 
Hoosic  River  flood-control  project  provides  flood  protection,  but  results  in 
harsh  conditions  for  fish  and  wildlife.  However,  a  restoration  plan  has 
been  proposed.  The  proposed  plan.  Section  1135  Ecosystem  Restoration,  is 
intended  to  restore  fisheries  and  wildlife  habitat  to  the  Hoosic  River  proj¬ 
ect.  Restoration  ideas  include  reshaping  of  the  trapezoidal  riprap  reach 
and  modifying  drop  structures  so  that  they  are  more  conducive  to  fish 
passage.  The  idea  is  to  add  a  winding  low-flow  channel  within  the  trape¬ 
zoidal  section  to  create  habitat.  Pools  and  riffles  would  be  formed  with  the 
low-flow  channel. 

The  proposed  modifications  investigated  in  this  report  were  changes  to  the 
downstream  5,000  ft  of  modeled  reach.  The  proposed  restoration  modifi¬ 
cations  begin  at  the  end  sill  of  drop  structure  No.  2  (sta  50+00)  and  ends 
immediately  upstream  of  the  Lime  Street  Bridge  (sta  1+50).  Numerical 
flow  models  were  used  to  evaluate  the  hydraulic  performance  of  two 
restoration  alternatives.  The  flow  velocities  and  water-surface  elevations 
resulting  from  the  proposed  modifications  were  computed  and  compared 
with  those  that  the  as-built  design  produced. 

A  temperature  model  was  used  by  the  NAN  to  size  the  low-flow  channel. 
These  results  indicated  that  a  16-ft-wide  low-flow  channel  should  be  incor¬ 
porated  into  the  restoration  plans.  Therefore,  a  16-ft  wide  by  2-  to  2.5-ft- 
deep  low-flow  channel  was  added  to  the  “Water  Resources  Development 
Project  Section  1135  Ecosystem  Restoration,  Stream  Temperature  Alter¬ 
native”  designs  provided  by  the  NAN.  The  low-flow  channel,  which  had  iV 
on  3H  side  slopes,  was  cut  into  the  existing  bed.  The  low-flow  channel 
meandered  from  the  right-  to  the  left-side-slope  toe.  The  existing  condi¬ 
tions  model,  discussed  in  the  previous  chapter,  was  modified  to  construct 
models  of  the  two  proposed  plans  titled  Alternative  A  and  Alternative  B. 
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Alternative  A  model 

Geometry  and  computational  mesh 

The  Alternative  A  restoration  design  had  a  meandering  low-flow  channel 
coupled  with  riffles  and  pools  at  the  crossings  and  bends.  Alternative  A 
called  for  replacing  the  existing  drop  structures  with  ramps  created  by 
Ailing  in  the  drops  with  stones.  These  ramps  had  an  incised  channel  which 
would  contain  the  daily  low  flows.  These  modifications  to  the  existing  drop 
structures  were  minimal.  The  existing  concrete  would  be  covered  with 
stacked  stones  forming  a  sloped  transition  across  the  drop  structure.  Each 
rock  ramp  was  240  ft  long  and  placed  on  a  3-percent  slope.  These  rock 
ramps  were  used  to  fill  in  drop  structures  No.  3  (crest  at  sta  36+54)  and 
No.  4  (crest  at  sta  20+67).  The  typical  modification  is  shown  in  Figure  56. 

The  computational  mesh  used  to  simulate  the  design  discharge  in  the 
Alternative  A  channel  is  shown  in  Figure  57.  The  mesh  had  10,496  nodes 
and  18,807  elements.  The  number  of  nodes  and  elements  was  less  than 
that  used  to  describe  the  as-built  or  existing  conditions  models  because  the 
rock  ramps  did  not  require  as  much  resolution  as  the  drop  structures  of 


Figure  56.  Typical  drop  structure  modification  used  in  the  restoration  alternatives. 
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Figure  57.  Alternative  A,  computationai  mesh  detaiis. 
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the  previous  models.  The  Alternative  A  model  included  the  concrete  high- 
velocity  channel  upstream  of  the  modified  reach  because  the  boundary 
conditions  were  established  and  the  finite  element  mesh  was  already 
constructed.  It  was  desirable  to  have  models  of  the  entire  reach  for  each 
configuration;  the  as-built,  existing,  Alternative  A,  and  Alternative  B,  in 
case  future  studies  investigate  discharges  and/or  conditions  other  than  the 
ones  described  in  this  report. 

Element  sizes  varied  throughout  the  domain,  but  the  primary  reaches  of 
the  trapezoidal  channel  were  described  with  elements  that  were  about 
20  ft  long  by  about  8  ft  wide.  The  width  of  the  low-flow  channel  was  usu¬ 
ally  described  using  2  elements.  The  side  slopes  of  the  low-flow  channel 
were  described  using  one  element  on  each  side  of  the  channel.  The  channel 
side  slopes  were  described  using  similar  sized  elements.  The  effects  of  the 
cobble  lateral  bars  and  riffles  were  modeled  with  the  elements’  roughness. 

The  modeling  process  consisted  of  time  marching  from  initial  conditions 
of  a  level  water-surface  and  quiescent  flow  to  steady  state  produced  by  the 
boundary  conditions.  The  mesh  was  converged  by  advancing  in  time  from 
the  original  mesh’s  steady-state  solution  to  a  steady-state  solution  with 
two  levels  and  three  levels  of  adaption.  The  steady-state  solution  with  two 
levels  of  adaption  did  not  significantly  differ  from  the  three  levels  and 
therefore  mesh  convergence  was  obtained  with  two  levels  of  adaption.  The 
two  levels  of  adaption  were  employed  on  the  low-flow  channel,  the  side 
slopes,  the  riffles,  and  the  rock  chutes  at  the  drop  structures.  The  auto¬ 
matic  adaption  process  produced  a  mesh  that  had  10,904  nodes.  This  is 
408  nodes  more  than  the  original  mesh. 

Boundary  conditions  and  modei  parameters 

The  inflow  discharges  and  tailwater  elevation  were  those  associated  with 
the  design  discharge  of  6,773  cfs  (100-year  frequency).  The  distribution  of 
inflows  is  provided  in  Table  1.  These  boundary  conditions  were  identical  to 
those  used  for  the  as-built  and  existing  conditions  models  since  the  pri¬ 
mary  purpose  of  this  study  was  to  evaluate  the  affects  that  restoration 
modifications  might  have  on  the  channel’s  ability  to  pass  the  design  dis¬ 
charge.  Simulations  were  run  to  evaluate  design  flow  conditions  resulting 
from  various  modifications  to  the  lower  end  of  the  project.  The  model  was 
used  to  determine  what  changes  the  restoration  alterations  would  make  to 
the  water-surface  elevations  and  flow  distributions  when  the  channel 
passed  the  design  discharge. 


ERDC/CHLTR-10-1 


69 


Roughness  coefficients  were  assigned  as  with  the  previous  models.  A  Man¬ 
ning’s  roughness  coefficient,  n,  of  0.014  was  used  for  concrete,  0.025  for 
the  existing  channel’s  gravel  bed,  0.020  for  the  assumed  sand  bottomed 
low-flow  channel,  and  0.030  for  the  sloped  grassy  banks.  The  cobbles  of 
the  riffles  and  the  stones  of  the  fill  slopes  at  the  drop  structures  were 
assigned  an  n  value  of  0.035. 

The  flow  at  the  model’s  outflow  boundary  was  subcritical  and  therefore  a 
tailwater  boundary  condition  was  set.  A  tailwater  elevation  at  the  outflow 
boundary  was  set  to  el  743.0.  This  is  the  water-surface  elevation  at  the 
normal  depth  produced  by  the  trapezoidal  cross  section  with  the  design 
discharge  of  6,773  cfs. 

Alternative  A  results 

Flow  velocities  at  the  two  ramps  that  are  currently  drop  structures, 
reached  values  near  19  fps.  These  ramps  will  require  large  size  stone  to 
resist  erosion  due  to  these  high  velocities.  There  was  very  little  difference 
in  the  depth-averaged  velocity  at  the  low-flow  channel  as  compared  with 
the  main  channel.  Although  the  velocity  was  slightly  slower  in  the  shallow 
area  near  the  waterline,  the  flow  distribution  at  the  design  discharge  was 
generally  uniform  across  a  cross  section.  Velocity  magnitude  contours  are 
shown  in  Figures  58-62.  Along  the  length  of  the  Alternative  A,  the  flow 
was  no  less  than  5  ft  deep  and  so  the  2-ft  added  depth  in  the  low-flow 
channel  did  not  change  the  depth-averaged  velocity  distribution 
significantly. 

Between  the  Tophet  Brook  confluence  and  drop  structure  No.  3  (crest  at 
sta  36+54),  the  depth-averaged  velocity  varied  from  about  9  fps  to  12.5  fps 
(Figures  58  and  59).  The  velocities  through  the  reach  from  drop  structure 
No.  3  to  drop  structure  No.  4  varied  from  about  6  fps  at  sta  31+50  to  16  fps 
near  sta  24+00  (Figures  59  and  60).  Downstream  of  drop  structure  No.  4 
the  largest  velocity  is  about  12.5  fps  near  sta  3+00  (Figures  60-62). 

Alternative  A  produced  a  water-surface  elevation  along  the  reach  between 
the  Tophet  Brook  confluence  and  drop  structure  No.  3  that  was  about  1  ft 
lower  than  that  produced  with  the  as-built  design  (Figures  63-66).  How¬ 
ever,  downstream  of  sta  36+00,  the  profiles  in  Figures  65-70  show  that 
Alternative  A  resulted  in  water-surface  elevations  that  were  as  large  as 
2.5  ft  higher  than  the  as-built  conditions. 
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Figure  58.  Alternative  A,  velocity  magnitude  contours,  sta  53+50  to  41+94.07. 
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Figure  59.  Alternative  A,  velocity  magnitude  contours,  sta  41+17  to  29+11. 


ERDC/CHLTR-10-1 


72 


ERDC/CHLTR-10-1 


73 


ERDC/CHLTR-10-1 


74 


55+00  54+00  53+00  52+00  51+00  50+00  49+00  48+00  47+00  46+00  45+00  44+00  43+00  42+00 

STATIONS  IN  FEET  ALONG  CENTER  LINE 


LEGEND 

— ALTERNATIVE  A  BED 
— AS-BUILT  BED 

— ALTERNATIVE  A  WATER  SURFACE 
— AS-BUILT  WATER  SURFACE 


DISCHARGE  STATIONS 
5229  CFS  103-150  TO  46-150 
6773  CFS  46-150  TO -1 -CIO 


PLAN  AND  LEFT  SIDE  PROFILES 

ALTERNATIVE  A 


NOTE:  SIDES  OF  CHANNEL  ARE  REFERENCED  TO  LOOKING  DOWNSTREAM 


STA  55+25.00  TO  STA  41+94.07 


Figure  63.  Alternative  A,  plan  and  profiles,  sta  55+25.00  to  41+94.07,  left  wall. 
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Figure  64.  Alternative  A,  pian  and  profiies,  sta  55+25.00  to  41+94.07,  right  waii. 
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Figure  65.  Alternative  A,  plan  and  profiles,  sta  41+94.07  to  27+19.78,  left  wall. 
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Figure  66.  Alternative  A,  pian  and  profiies,  sta  41+94.07  to  27+19.78,  right  waii. 
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Figure  67.  Alternative  A,  plan  and  profiles,  sta  27+19.78  to  12+14.29,  left  wall. 
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Figure  68.  Alternative  A,  plan  and  profiles,  sta  27+19.78  to  12+14.29,  right  wall. 
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Figure  69.  Alternative  A,  plan  and  profiles,  sta  12+14.29  to  0+00,  left  wall. 
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Figure  70.  Alternative  A,  pian  and  profiies,  sta  12+14.29  to  0+00,  right  waii. 
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Alternative  B  model 

Geometry  and  computational  mesh 

The  Alternative  B  called  for  cutting  the  concrete  of  the  existing  drop  struc¬ 
tures.  The  invert  at  the  crest  of  the  existing  drop  structures  would  be  cut 
4  ft.  The  remaining  drop  would  then  be  filled  with  rock  to  form  a  240-ft- 
long  ramp  having  a  1.5-percent  slope.  This  alternative  would  change  the 
entire  longitudinal  slope  of  the  Hoosic  River  between  drop  structure  No.  2 
(crest  at  sta  51+50)  and  No.  4  (crest  at  sta  20+67). 

The  Alternative  B  design  replaced  the  rectangular  cross  section  between 
sta  51+01.5  and  sta  44+50  with  a  trapezoidal  cross  section.  The  existing 
vertical  concrete  walls  through  this  reach  were  replaced  with  iV  on  3H 
side  slopes  armored  with  rock.  The  toe  of  the  side  slopes  were  at  the  base 
of  the  existing  walls  so  that  the  Alternative  B  channel  was  significantly 
wider  than  the  existing  channel  through  this  reach. 

A  low-flow  channel  similar  to  Alternative  A  was  cut  into  the  existing  chan¬ 
nel  bed  in  a  pattern  similar  to  Alternative  A.  Pools  and  riffles  were  formed 
as  the  low-flow  channel  crossed  from  side  to  side  of  the  bed.  The  primary 
difference  in  the  beds  of  Alternative  A  and  B  is  that  Alternative  B  removed 
the  drop  structures  and  that  the  resulting  bed  slope  was  steeper  than  the 
slopes  between  the  drop  structures  of  Alternative  A.  The  purpose  of  the 
existing  drop  structures  was  to  reduce  the  bed  slope  to  ensure  that  the  flow 
is  subcritical.  Alternative  B  lowered  the  existing  bed  along  the  upstream 
3,100  ft  (sta  49+00  to  sta  18+00)  of  the  restoration  project.  Alternative  B 
was  similar  to  Alternative  A  downstream  of  sta  18+00  where  the  low-flow 
channel  is  cut  into  the  existing  trapezoidal  channel. 

Alternative  B  called  for  “bank  plantings”  along  both  side  slopes.  The  resto¬ 
ration  line  drawings  showed  that  these  plantings  would  include  bushes 
and  trees.  The  bank  plantings  were  not  simulated  in  the  Alternative  B 
model.  The  U.S.  Army  Corps  of  Engineers’  guidance  requires  a  vegetation- 
free  zone  (Headquarters,  U.S.  Army  Corps  of  Engineers  (USACE)  2000) 
such  that  for  trapezoidal  channels,  trees  can  be  planted  no  closer  than  15  ft 
from  the  toe  of  the  outside  levee  (USACE  2009).  Therefore,  the  roughness 
coefficients  did  not  reflect  the  presence  of  these  plantings,  but  rather  were 
of  the  same  value  as  those  used  to  represent  the  roughness  of  Alternative  A 
side  slopes. 
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The  Alternative  B  mesh  had  10,962  nodes  and  19,775  elements.  Element 
sizes  and  mesh  resolution  were  similar  to  those  used  to  model  the  Alterna¬ 
tive  A  conditions.  Details  of  the  Alternative  B  mesh  are  shown  in  Figure  71. 
Two  levels  of  adaption  were  employed  on  the  low-flow  channel,  the  side 
slopes,  the  riffles,  and  the  rock  chutes  at  the  drop  structures.  Adaption 
resulted  in  an  addition  of  347  nodes  to  the  original  mesh  to  establish  the 
steady-state  solution. 

Boundary  conditions  and  modei  parameters 

The  boundary  conditions  were  identical  to  those  used  to  model  Alterna¬ 
tive  A.  The  design  discharge  of  6,773  cfs  was  the  cumulative  inflow  and  a 
tailwater  of  el  743.0  was  set  at  the  outflow. 

The  roughness  coefficients  were  the  same  as  those  used  in  the  Alternative 
A  model.  A  Manning’s  roughness  coefficient,  n,  of  0.014  was  used  for  con¬ 
crete,  0.025  for  existing  channel’s  gravel  bed,  0.020  for  the  assumed  sand 
bottomed  low-flow  channel,  0.030  for  the  sloped  grassy  banks,  and  0.035 
for  the  riffles  and  the  stones  of  the  fill  slopes  at  the  drop  structures. 

Aiternative  B  resuits 

Contours  of  velocity  magnitudes  with  Alternative  B  are  shown  in  Fig¬ 
ures  72-76.  The  velocity  immediately  downstream  of  the  Tophet  Brook 
confluence  was  not  significantly  different  between  the  two  alternatives. 
However,  the  Alternative  B  was  wider  between  drop  structure  No.  2  and 
sta  44+50  and  the  Tophet  Brook  flow  crossed  the  main  channel  flow  more 
than  the  Alternative  A  flow.  The  wider  cross  section  of  the  Alternative  B 
design  resulted  in  depth-averaged  velocities  being  about  0.5  fps  lower  than 
those  produced  with  Alternative  A. 

The  Alternative  B  had  milder  longitudinal  slopes  (1.5  percent)  over  the 
existing  drop  structures  than  Alternative  A  (3  percent)  and  so  the  veloci¬ 
ties  were  slightly  lower.  The  maximum  velocity  at  the  drop  structures  was 
about  16  fps  with  Alternative  B  compared  to  19  fps  with  the  Alternative  A. 
The  Alternative  B  flow  at  the  existing  drop  structures  is  about  0.5  ft  to 
1.0  ft  deeper  than  Alternative  A. 
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Details  at  drop  structure  4 
and  stilling  basin 


Details  at  drop  structure  3  and 
stilling  basin 


Details  at  drop  structure  2  and 
stilling  basin 


Figure  71.  Alternative  B,  computational  mesh  details. 
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Figure  72.  Alternative  B,  velocity  magnitude  contours,  sta  53+50  to  41+94.07. 
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The  overall  longitudinal  bed  slope  of  Alternative  B  was  steep  enough  to 
produce  flow  that  is  near  the  critical  condition.  The  critical  flow  condition 
is  where  the  Fronde  number  nears  unity.  Channels  flowing  near  the  critical 
condition  will  experience  unstable  flow  with  an  undulating  water  surface. 
The  numerical  model  results  are  steady-state  conditions  and  do  not  reflect 
the  instability  associated  with  flow  near  critical  (Fronde  number  near 
unity).  The  design  manual  EM  1110-2-1601  (USAGE  1991)  notes  that  when 
flow  nears  the  critical  condition  “a  relatively  large  change  of  depth  may 
occur  with  a  very  small  variation  of  specific  energy.  Elow  in  this  region  is 
unstable  and  excessive  wave  action  or  undulations  of  the  water  surface 
may  occur.” 

Corps  of  Engineers’  guidance  (USAGE  1991)  recommends  that  channels  be 
designed  to  avoid  Eroude  numbers  between  0.86  and  1.13.  It  is  noted  that 
the  numerical  model  results  presented  in  this  report  do  not  consider  the 
undulations  of  the  water  surface  that  may  occur. 

Alternative  B  produced  water-surface  elevations  lower  than  the  as-built 
design  because  the  channel  was  bed  was  lower  and  steeper  (and  wider  in 
the  reach  where  side  slopes  replaced  vertical  walls).  The  profiles  provided 
in  Eigures  77-84  show  that  the  water-surface  elevations  from  drop  struc¬ 
ture  No.  2  (crest  at  sta  51+00)  down  to  drop  structure  No.  4  (crest  at  sta 
20+67)  were  significantly  lower  than  the  as-built  design.  This  is  because 
the  bed  of  the  Alternative  B  was  lower  than  the  as-built  design.  Down¬ 
stream  of  sta  17+00,  the  water-surface  profiles  with  the  Alternative  B  were 
essentially  the  same  as  those  produced  with  Alternative  A  since  the 
designs  were  similar  through  this  reach. 
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Figure  77.  Alternative  B,  plan  and  profiles,  sta  55+25.00  to  41+94.07,  left  wall. 
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Figure  78.  Alternative  B,  plan  and  profiles,  sta  55+25.00  to  41+94.07,  right  wall. 
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Figure  79.  Alternative  B,  plan  and  profiles,  sta  41+94.07  to  27+19.78,  left  wall. 
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Figure  80.  Alternative  B,  plan  and  profiles,  sta  41+94.07  to  27+19.78,  right  wall. 
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Figure  81.  Alternative  B,  plan  and  profiles,  sta  27+19.78  to  12+14.29,  left  wall. 
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Figure  82.  Alternative  B,  plan  and  profiles,  sta  27+19.78  to  12+14.29,  right  wall. 
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Figure  83.  Alternative  B,  plan  and  profiles,  sta  12+14.29  to  0+00,  left  wall. 
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Figure  84.  Alternative  B,  pian  and  profiies,  sta  12+14.29  to  0+00,  right  waii. 
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6  Summary 

This  study  determined  the  design-flow  conditions  in  the  man-made  reach 
of  the  Hoosic  River  that  forms  the  flood-control  project  for  Adams,  MA. 
Flow  conditions  were  simulated  with  the  2-D  depth-averaged  flow  code, 
ADH.  These  simulations  were  conducted  to  determine  water-surface  ele¬ 
vations  and  flow  velocities  along  the  river  from  upstream  of  the  weir  near 
sta  205+96.03  where  the  project  is  a  rectangular  concrete  channel  to  the 
Lime  Street  Bridge  where  it  is  a  trapezoidal  riprap  channel. 

The  model  was  validated  by  comparing  with  previously  published  physical 
model  data  for  the  Hoosic  River  Flood-Control  Channel  through  North 
Adams,  MA.  Simulation  results  showed  that  the  2-D  ADH  model  is  an 
appropriate  tool  to  use  in  evaluating  flow  conditions  in  the  Hoosic  River 
Flood-Control  Channel. 

As-built  and  existing  conditions  were  modeled  primarily  to  determine  the 
consequences  of  the  sediment  deposition  in  the  lower  reach.  Water-surface 
elevations  produced  during  the  100-year  frequency  flow  were  compared. 
The  water-surface  elevations  produced  with  the  existing  conditions  were 
more  than  2  ft  in  areas  of  the  channel  where  deposits  of  sediment 
obstructed  the  flow. 

The  primary  purpose  of  the  study  was  to  evaluate  proposed  channel  rest¬ 
oration  modifications  to  the  lower  end.  The  model  determined  whether 
restoration  alterations  changed  water-surface  elevations  associated  with 
the  design  discharge.  Restoration  designs  included  a  16-ft-wide  by  2-  to 
2.5-ft-deep  low-flow  channel  with  pools  and  riffles.  Two  restoration  con¬ 
figurations  were  studied. 

Alternative  A,  downstream  of  sta  36+00,  produced  water-surface  eleva¬ 
tions  that  were  in  some  areas  2.5  ft  higher  than  the  as-built  conditions  at 
the  design  discharge.  However,  some  reaches  of  Alternative  A  resulted  in 
water-surface  elevations  that  were  lower  than  the  as-built  conditions. 

Although  Alternative  B  produced  water-surface  elevations  that  were  lower 
than  the  as-built  conditions,  the  design  discharge  produced  flow  near  criti¬ 
cal  depth.  This  opposes  Corps  of  Engineers’  guidance  (USAGE  1991)  which 
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recommends  that  channels  be  designed  to  avoid  Fronde  numbers  near 
unity.  The  numerical  model  results  presented  in  this  report  do  not  con¬ 
sider  the  undulations  of  the  water  surface  that  may  occur.  Alternative  A 
provided  better  hydraulic  conditions  than  Alternative  B  because  the  later 
resulted  in  flow  near  critical  conditions. 
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